Many-body large polaron optical conductivity in SrTii_ x Nb z 03 
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Recent experimental data on the optical conductivity of niobium doped SrTiOs are interpreted 
in terms of a gas of large polarons with effective coupling constant a e ff « 2. The theoretical 
approach takes into account many-body effects, the electron-phonon interaction with multiple LO- 
phonon branches, and the degeneracy and the anisotropy of the Ti t2 9 conduction band. Based 
on the Frohlich interaction, the many-body large-polaron theory provides an interpretation for the 
essential characteristics, except - interestingly - for the unexpectedly large intensity of a peak at 
~ 130 meV, of the observed optical conductivity spectra of SrTii-^NbxOa without any adjustment 
of material parameters. 



I. INTRODUCTION 

The infrared optical absorption of perovskite-type ma- 
terials, in particular, of copper oxide based high-T c super- 
conductors and of the manganites has been the subject 
of intensive investigations [3TfT2], Insulating SrTiC>3 has 
a perovskite structure and manifests a metal-insulator 
transition at room temperature around a doping of 
0.002% La or Nb per unit cell [13]. At low doping con- 
centrations, between 0.003% and 3%, strontium titanate 
reveals a superconducting phase transition [T3] below 
0.7 K. Various optical experiments [T3J IT5HTU] show a 
mid-infrared band in the normal state optical conduc- 
tivity of doped SrTiC>3 which is often explained by po- 
laronic behavior. In the recently observed optical con- 
ductivity spectra of Ref. [T7j, shown in Fig. 1, there is 
a broad mid-infrared optical conductivity band starting 
at a photon energy of Ml ~ 100 meV, which is within 
the range of the LO-phonon energies of SrTii_ x Nb x C>3. 
The peaks/shoulders of the experimental optical conduc- 
tivity band at Kl ~ 200 to 400 meV resemble the peaks 
provided by the mixed plasmon-phonon excitations as de- 
scribed in Ref. [20 . Based on the experimental data, the 
authors deduce a coupling constant 3 < a < 4 and con- 
clude the mid-infrared peaks to originate from large po- 
laron formation. The high and narrow peaks positioned 
at the lower frequencies with respect to the mid-infrared 
band are attributed in Ref. [17] to the optical absorption 
of the TO-phonons. 

There are different types of polaron states in solids. In 
the effective mass approximation for the electron placed 
in a continuum polarizable medium, a so-called large or 
continuum polaron can exist. Large polaron wave func- 
tions and the corresponding lattice distortions spread 
over many lattice sites. Due to the finite phonon frequen- 
cies the ion polarizations can follow the polaron motion 
if the motion is sufficiently slow. Hence, large polarons 
with a low kinetic energy propagate through the lattice 
as free electrons but with an enhanced effective mass. 
When the polaron binding energy is larger than the half 
bandwidth of the electron band, all states in the Bloch 
bands are 'dressed' by phonons. In this strong-coupling 
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FIG. 1: Optical conductivity of SrTii-^Nb^Oa for 0.1% (grey 
curves), 0.2% (blue curves), 0.9% (green curves) and 2% (pink 
curves) at 300 K (panel a) and 7 K (panel b). For clarity, 
the mid-infrared conductivities of x = 0.1% and 0.2% are 
magnified by a factor 4. (From Ref. [17].) 



regime, the finite electron bandwidth becomes important, 
so the continuum approximation cannot be applied. In 
this case the carriers are described as "small" or discrete 
(lattice) polarons that can hop between different states 
localized at lattice sites. A key distinction between large 
and small polarons is then the radius of the polaron state. 
For large polarons, that radius substantially exceeds the 
lattice constant, while for small polarons it is compara- 
ble to the lattice constant. A review of the properties 
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of large and small polarons can be found, e. g., in Refs. 
[2TJ [22]. In the theory of "mixed" polarons [25H26] the 
states of the electron-phonon system are composed of a 
mixture of large and small polaron states. 

Polaron states are formed due to the electron-phonon 
interaction, which is different in the cases of large and 
small polarons. For a large polaron, the electron-phonon 
interaction is provided by a macroscopic (continuum) po- 
larization of the lattice. This interaction is characterized 
by the coupling constant a introduced by Frohlich [27], 
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where and £q are, respectively, the high-frequency and 
the static dielectric constants, e is the electron charge, 
ljl is the longitudinal optical (LO) phonon frequency 
in the Brillouin zone centre, and is the band elec- 
tron (or hole) mass. The large-polaron coupling con- 
stant is thus expressed through macroscopic observable 
parameters of a polarizable medium. On the contrary, 
the electron-phonon interaction for a small polaron is de- 
scribed through microscopic parameters. 

The nature of the polaron states in SrTii-^NbxOs is 
not yet clear. Previous optical measurements on stron- 
tium titanate were interpreted in terms of small polarons 
[21HI28] . However, that assumption contradicts the inter- 
pretation of transport measurements [29) . which rather 
support the large-polaron picture. Also the heat capac- 
ity measurements |30j , provide effective masses similar to 
those of large polarons. In Ref. [15], the experimental re- 
sults of Ref. |18j on the temperature-dependent plasma 
frequencies in SrTii-^Nb^Oa were interpreted within the 
theory of mixed polarons [23-26 . Thermoelectric power 
measurements [29 i have shown that the density-of-states 
masses increase with increasing temperature, which can 
be explained by a theory of mixed polarons |23j . It has 
been supposed [16] that the polaron optical conductiv- 
ity in SrTii_ ;E Nb ;E 03 is probably provided by mixed po- 
larons. A possible coexistence of large and small mass 
polarons has been suggested in Ref. [31. In Ref. [32] . 
coexistence of small and large polarons in the same solid 
is invoked to interpret experimental data on the optical 
absorption in oxides. 

The key question is to determine the type of polarons 
that provide the mechanism of the polaron optical con- 
ductivity in SrTii_j;Nb K 03. The optical response of large 
polarons in various approximations was studied, e. g., 
in Refs. [33-36]. The same problem for the small po- 
laron was investigated in [2E1 (S3 • ^ n the large-polaron 
theory, the optical absorption is provided by transitions 
(with 0, 1, . . . phonon emission) between different contin- 
uum electron states. In the small-polaron theory, the 
optical absorption occurs when the self-trapped carrier 
is induced to transfer from its localized state to a local- 
ized state at an adjacent site, with emission of phonons. 
Because of the different physical mechanisms involved, 
the optical conductivity spectra of large and small po- 
larons are different from each other. In the large-polaron 



theory the polaron optical conductivity behaves at high 
frequencies Q as a power function (exfi -5 / 2 ). In the 
small-polaron theory, the polaron optical conductivity at 
high frequencies decreases much faster than for large po- 
larons: as a Gaussian exponent. Therefore the analysis of 
optical measurements can shed some light on the afore- 
said question on the type of polarons responsible for the 
optical conductivity in SrTii-^Nb^Oa. 

The polaron optical conductivity band of 
SrTi 1 _ :c Nb :E 03 occupies the mid-infrared range of 
the photon energies Ml < 1 eV, and the threshold for 
interband electron-hole transitions lies at the band gap 
energy, which is around 3.3 eV in SrTii_ x Nba;03 [17]. 
Therefore interband transitions do not interfere with 
the polaron optical conductivity. Other mechanisms 
of electron intraband scattering (for example, electron- 
phonon interaction with acoustic phonons and/or 
electron or hole transitions from impurity centers) may 
be manifested together with the polaron mechanism in 
the energy range f&l < 1 eV. The treatment of those 
mechanisms is, however, beyond the scope of the present 
investigation. 

We can make some preliminary suggestions concern- 
ing the dominating mechanism of the mid-infrared op- 
tical conductivity in the Nb doped strontium titanate. 
The low-frequency edge of the mid-infrared band in 
SrTii-aNbaOs at a low temperature (T = 7 K) lies in the 
range of the LO- phonon energies obtained in [18] . The 
maximum of the mid-infrared band lies relatively close to 
this low-frequency edge (the difference in frequency be- 
tween the low-frequency edge and the maximum of the 
mid-infrared band is comparable to the LO-phonon fre- 
quencies in SrTii-^NbxOa). This behavior is charac- 
teristic of large-polaron optical conductivity rather than 
of small-polaron optical conductivity. Indeed, the max- 
imum of the small polaron optical conductivity band is 
expected to be shifted to considerably higher frequencies 
with respect to the low-frequency edge of the polaron 
optical conductivity band (see, e.g., Ref. [37]). Also, 
at sufficiently high frequencies, the experimental mid- 
infrared band from Ref. [T7] decreases with increasing 
n rather slowly, which is characteristic for large-polaron 
optical conductivity rather than for small-polaron op- 
tical conductivity. We therefore can suggest that the 
large-polaron picture is the most appropriate for the in- 
terpretation of the mid-infrared band of SrTii-^Nb^Os 
observed in Ref. |17j . 

In order to interpret the mid-infrared band of 
the experimental optical conductivity spectra of 
SrTii_a;Nb x 03 [T7] in terms of polarons, we calcu- 
late the large-polaron optical conductivity spectra for 
SrTii_a;Nb K 03 using the model for the optical conductiv- 
ity of a large-polaron gas developed in Ref. [20] , adapted 
to take into account multiple LO-phonon branches |38) . 
The degeneracy and the anisotropy of the conduction 
band in SrTii-^Nb^Oa are taken into account. 

The paper is organized as follows. In Sec. [IT] we de- 
scribe the theoretical formalism for the polaron optical 
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conductivity. In Sec. Ill the calculated polaron optical 



conductivity spectra are discussed and compared with 



the experiment [T7]. Sec. IV contains conclusions. 



operators for a phonon of the j-th branch with the mo- 
mentum q, Uij- are the LO-phonon frequencies (approx- 
imated here as non-dispersive), and V is the volume of 
the crystal. The polaron interaction amplitude Vqj is 



II. OPTICAL CONDUCTIVITY OF A GAS OF 
LARGE POLARONS 

The optical absorption spectra of SrTii-^Nb^Oa are 
sensitive to the doping level [T7]. Therefore a many- 
polaron description is in order. In our context, "many- 
polaron description" means an account of many-electron 
effects on the optical conductivity of a polaron gas. 
These effects include the influence of the electron-electron 
Coulomb interaction (which leads to screening effects) 
and of the Fermi statistics of the polaron gas on the 
optical conductivity spectra. In the low-density limit, 
those many-body effects are not important, and the op- 
tical conductivity of a polaron gas is well described by the 
optical conductivity of a single polaron multiplied by the 
electron density. The scope of the present study embraces 
a wide range of electron densities for which the single- 
polaron approach is, in general, insufficient. As shown 
below, even at the lowest electron density involved in the 
experiment |17) . the shape and magnitude of the optical 
conductivity spectrum is strongly affected by many-body 
effects. 

We wish to compare the experiments of Ref. [T7], in 
particular the observed mid-infrared band, to the the- 
oretical optical conductivity of a gas of large polarons. 
For that purpose we use the many-body large polaron 
approach of Refs. [501 13H] > which takes into account the 
electron-electron interaction and the Fermi statistics of 
polarons. 

Refs. [20., ,38] are limited to the study of weak-coupling 
polarons. Up to a ~ 3, which includes the case of 
SrTii-^NbzOa, the weak coupling approximation can 
be expected to describe the main characteristics of the 
many-polaron optical response (see, e.g., Refs. [20H22] ). 
In Ref. [35] a generalization of Ref. [5D] is presented 
that takes into account the electron-phonon interaction 
with multiple LO-phonon branches as they exist, e. g., 
in complex oxides. For a single polaron, effects related to 
multiple LO-phonon branches were investigated in Ref. 
[4Ti] . The starting point for the treatment of a many- 
polaron system is the Frohlich Hamiltonian 
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where ctj is a dimensionless partial coupling constant 
characterizing the interaction between an electron and 
the j-th LO-phonon branch. The electron-electron inter- 
action is described by the Coulomb potential energy 
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Optical phonons in SrTiOa show a considerable disper- 
sion (see, e. g., Ref. [H] and references therein). The ef- 
fect of the phonon dispersion is a broadening of features 
of the polaron optical conductivity band. The magni- 
tude of the broadening is characterized by the dispersion 
parameter Aw of the optical phonons, that contribute 
to the integrals over q entering the polaron optical con- 
ductivity. In a polar crystal with a single LO-phonon 
branch, that range of convergence is approximately qo = 

1/2 

{mbLOLo/K) . For SrTi03, taking w;,o = rnaxjwLj}, 
we obtain qo ss 1.02 x 10 9 m _1 . The boundary of the 
Brillouin zone tt/oq in SrTiOa (where the lattice con- 
stant ao R3 0.3905 nm) is at 8 x 10 9 m _1 . Therefore the 
integration domain for the relevant integrals is one order 
smaller than the size of the Brillouin zone. In the region 
< q < qo, the dispersion parameter of the LO-phonon 
frequencies, Aw, is a few percent of lolj- Consequently, 
Aw is very small compared with the characteristic width 
of the polaron band. Therefore, in the present treatment, 
we apply the approximation of non-dispersive phonons. 

For a description of a polarizable medium with n 
optical-phonon branches, we use the model dielectric 
function HI E3J 
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where cJ CT (ck.o-) are the creation (annihilation) opera- 
tors for an electron with momentum k and with the spin 
z-projection er, alj (a q ,j) are the creation (annihilation) 



whose zeros (poles) correspond to the LO(TO) phonon 
frequencies u>L,j if^T j ) ■ This dielectric function is the re- 
sult of the straightforward extension of the Born-Huang 
approach |H] to the case where more than one optical- 
phonon branch exists in a polar crystal. The Born-Huang 
approach and its extension |42j generate expressions for 
the macroscopic polarization induced by the polar vi- 
brations, and for the corresponding electrostatic poten- 
tial. This electrostatic potential is a basis element of the 
Hamiltonian of the electron-phonon interaction. In Ref. 
[42] , the Hamiltonian of the electron-phonon interaction 
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has been explicitly derived with the amplitudes 
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Using Eqs. (|3| and ([6]) with the dielectric function 
we arrive at the following set of linear equations for the 
coupling constants otj (j = 1, . . . , n): 
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Knowledge of the band mass, of the electronic dielectric 
constant and of the LO- and TO-phonon frequencies 
is sufficient to determine the coupling constants (Xj tak- 
ing into account mixing between different optical-phonon 
branches. In the particular case of a single LO-phonon 
branch, Eq. ^ is reduced to ([I]). 

In order to describe the optical conductivity of a po- 
laron gas, we refer to the work [46] . where the Mori- 
Zwanzig projection operator technique has been used to 
rederive the path-integral result of Ref. [45] and the 
impedance of Ref. [33]. We repeat the derivations of 
Ref. 46 with the replacement of single-electron func- 
tions by their many-electron analogs. For example, e lq r 
in the Hamiltonian of the electron-phonon interaction is 
replaced by the Fourier component of the electron density 
for an TV-electron system, 
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As a result, we arrive at a formula which is structurally 
similar to the single-polaron optical conductivity [3"3"ll46| . 
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where n = N/V is the carrier density, and \ (O) is the 
memory function. The same many-electron derivation as 
in the present work, to the best of our knowledge, was 
first performed for the polaron gas in 2D in Ref. [47] in 
the weak electron-phonon coupling limit. 

In Refs. [33 EH] the single-polaron memory function 
was calculated starting from the all-coupling Feynman 
variational principle |48) . For a many-polaron system, an 
effective all-coupling extension of that variational princi- 
ple has not been worked out yet. In the present treat- 
ment, we restrict ourselves to the weak-coupling approx- 
imation for the electron-phonon interaction to derive the 
memory function. In this approximation, the memory 
function x (fi) is similar to that of Ref. [47 , with two 
distinctions: (1) the electron gas in the present treatment 
is three-dimensional, (2) several LO phonon branches are 
taken into account. The resulting form of the memory 



function is 
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where j3 = 1/ (ksT). The dynamical structure factor 
S (q, t) is proportional to the two-point correlation func- 
tion (cf. Ref. [20]), 
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To obtain x (^) to order a it is sufficient to perform 
the averaging in the correlation function (11) using the 
Hamiltonian ^ without the electron-phonon interaction 
and keeping the electron-electron interaction term U e - e . 
We calculate the dynamical structure factor (11) ex- 



tending the method [20] to nonzero temperatures. In Ref. 
[20) . the key advantage of the many-polaron variational 
approach [15] is exploited: the fact that the many-body 
effects are entirely contained in the dynamical structure 
factor S (q, t). The structure factor can be calculated us- 
ing various approximations. Terms of order of | Vq^- 1 2 are 
automatically taken into account in the memory function 
(10). Consequently, up to order a for <j(Q), it is suffi- 



cient to calculate S (q, t) without the electron-phonon 
coupling. In Ref. [5D], S (q, t) was calculated within two 
different approximations: (i) the Hartree-Fock approxi- 
mation, (ii) the random-phase approximation (RPA). As 
shown in Ref. [20 , the RPA dynamical structure fac- 
tor, contrary to the Hartree-Fock approximation, takes 
into account the effects both of the Fermi statistics and 
of the electron-electron interaction on the many-polaron 
optical-absorption spectra. 

The dynamical structure factor is expressed through 
the density-density Green's functions defined as 

poo 
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In terms of G (q, O) and G R (q, O), the memory function 
( 10 ) takes the form: 
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Taking into account the Coulomb electron-electron in- 
teraction within RPA, the retarded Green's function 
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G R (q, fi) is given by 



III. THEORY AND EXPERIMENT 



G*(q,fi) 



hvp^ (q, n) 



4-Tre 2 



(q,0) 



(15) 



where (q, f2) is the polarization function of the free 
electron gas, see, e.g., [50] 
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with the electron average occupation numbers /k.er- The 
function Q (q, J7) is obtained from G R (q, f2) using the 
exact analytical relation 



(1- 
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and the Kramers-Kronig dispersion relations for (q, Q). 

The above expressions are written for an isotropic 
conduction band. However, the conduction band of 
SrTii-zNbzOs is strongly anisotropic and triply degen- 
erate. The electrons are doped in three bands: d xy , d yz 
and d xz , which all have their minima at k = 0. Each 
of these bands has light masses along two direction (x 
and y for d xy , etc.) and a heavy mass along the third 
direction. While each electron has a strongly anisotropic 
mass, the electronic transport remains isotropic due to 
the fact that 2 light masses and 1 heavy mass contribute 
along each crystallographic axis. 

The anisotropy of the electronic effective mass of the 
conduction band can be approximately taken into ac- 
count in the following way. We use the averaged inverse 
band mass 



1 



3 \ m T 



1 

TO,, 



1 

m, 



and the density-of-states band mass 
m D = {m x m y m z ) 



(18) 



(19) 



The mass mo appears in the prefactor of the linear term 
of the specific heat. Comparing the mass Tod obtained 
from the experimental specific heat [301 EI] with the mass 
fhb obtained using optical spectral weights [T7] reveals 
the mass ratio of the heavy and light bands to be about 
27. The expression ( 18 ) replaces the bare mass m& in the 
optical conductivity and in the memory function ( 14 1. 
The polarization function of the free electron gas ( 16 ) is 



calculated with the density-of-states mass instead of 
mi,. The band degeneracy is taken into account through 
the degeneracy factor which is equal to 3, both in the 
polarization function and in the normalization equation 
for the chemical potential. The reduction of the polaron 
optical conductivity band due to screening with band de- 
generacy turns out to be less significant than without 
band degeneracy. 



A. Material parameters 

Several experimental parameters characterizing 
SrTii-aNbjcOs are necessary for the calculation of 
the large-polaron optical conductivity (see, e.g., Refs. 
[ISIHB]): the LO- and TO-phonon frequencies, the elec- 
tron band mass, and the electronic dielectric constant 

The electronic dielectric constant can be obtained us- 
ing reflectivity spectra of SrTii-^Nb^Oa. At T = 10 
K, the reflectivity of SrTii_ a .Nb x 3 is R « 0.16 for 

fa 5000 cm -1 . The electronic dielectric constant can 
be approximated using the expression 



R(n) 



v^P)-i 



(20) 



and assuming that fl = 5000 cm 1 is a sufficiently high 
frequency to characterize the electronic response. From 
pOb it follows that for SrTix^Nb^Oa, e Q 



5.44. 



In order to determine the optical-phonon frequencies, 
we use the experimental data from available sources [171 
[18] . In Ref. [17] . three infrared active phonon modes 
are observed at room temperature: at 11.0 meV, 21.8 
meV and 67.6 meV. With decreasing temperature, the 
lowest-frequency infrared-active phonon mode shows a 
strong red shift upon cooling, and saturates at about 2.3 
meV at 7 K. Those infrared-active phonon modes are 
associated with the polar TO-phonons. The TO-phonon 
frequencies determined in Ref. 18J for SrTii-^Nb^Oa 
with x = 0.9% at T = 300 K are 100 cm" 1 , 175 cm" 1 
and 550 cm -1 . The corresponding TO-phonon energies 
are 12.4 meV, 21.7 meV and 68.2 meV. 

Refs. [17] and [18] are used as sources for phonon pa- 
rameters. In Ref. [18], the TO-phonon frequencies are 
calculated on the basis of reflectivity measurements using 
a model dielectric function to fit experimental data. In 
Ref. [T7], the TO-phonon frequencies are obtained from 
an analysis of both reflectivity and transmission spectra, 
using inversion of the Fresnel equations of reflection and 
transmission coefficients and the Kramers-Kronig trans- 
formation of the reflectivity spectra. The TO-phonon en- 
ergies reported in Refs. [T7] and [TB] are in close agree- 
ment. This confirms the reliability of both experimen- 
tal data sources [T7J [TH] . The values of the TO-phonon 
frequencies used in our calculation are taken from the 
experiment [T7] because they are directly related to the 
samples of SrTii-^Nb^Oa for which the comparison of 
theory and experiment is made in the present work. 

The TO phonon frequencies from Ref. [17] can be used 
when they are complemented with corresponding LO 
phonon frequencies. However, Ref. [T7] does not contain 
data of the LO-phonon frequencies. In the present calcu- 
lation we use the LO phonon frequencies from Ref. [15] . 

The averaged band mass ( 18 1 is taken to be to/, = 
0.8lTO e (where m e is the electron mass in vacuum) ac- 
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TABLE I: Optical-phonon frequencies and partial coupling constants of doped strontium titanate 
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67.1 
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98.7 
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98.7 
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98.7 


98.7 


a 3 


1.582 


1.582 


1.582 


1.580 


1.576 


1.578 


1.570 


1.574 


a e ft 


2.06 


2.01 


2.06 


2.01 


2.05 


2.01 


2.03 


2.01 



cording to experimental data from Ref. [52] ■ Using the 
ratio of the heavy mass (m z ) to the light mass (m x — 
m y)i rn z /m x — 27, we find the density-of states band 
mass mjj rj 1.65m e . 

The TO- and LO- phonon frequencies and the resulting 
partial coupling constants calculated using the mass fhb 
are presented in Table 1. 

The effective coupling constant in Table 1 is deter- 
mined following Ref. [10], as a sum of partial coupling 
constants a,-, 



(21) 



The result a e ff ~ 2 shows that the electron-phonon cou- 
pling strength in SrTi 1 _ :c Nb ;! ,03 lies in the intermedi- 
ate to weak coupling range, and the conditions for small 
polaron formation are not fulfilled. This analysis indi- 
cates that the large-polaron picture - rather than the 
small-polaron description is suitable for the interpreta- 
tion of the mid-infrared band of the optical conductivity 
of SrTir^Nb.Oa. 

We use the actual electron densities for the samples 
studied in Ref. [T7] based on the unit cell volume (59.5 
cubic angstrom) and the chemical composition [x is the 
doping level). These carrier densities (see Table 2) are 
confirmed by measurements of the Hall constants. 



TABLE II: Electron densities of SrTii^Nb^Og 



x(%) 


n (cm -3 ) 


0.1 


1.7 x 10 19 


0.2 


3.4 x 10 19 


0.9 


1.5 x 10 20 


2.0 


3.4 x 10 20 



B. Optical conductivity spectra 

We calculate the large-polaron optical conductivity 
spectra for SrTii-^NbzOa using the approach of Ref. 
[213] as adapted in Ref. [35] to take into account multiple 
LO-phonon branches. We also include in the numerical 
calculation the TO-phonon contribution to the optical 
conductivity, described by an oscillatory-like model di- 
electric function (see, e.g., Ref. [15]): 



Re a TO (O) = — 



^2 , 2 



(22) 



where the weight coefficients <7qj and the damping pa- 
rameters 7j for each j-th TO-phonon branch are ex- 
tracted from the experimental optical conductivity spec- 
tra of Ref. [17]. The polaron-and the TO-phonon opti- 
cal responses are treated as independent of each other. 
Consequently the polaron-Q and TO-phonon (22) con- 
tributions enter the optical conductivity additivcly. 

Following the procedure described above using the ma- 
terial parameters discussed above, we obtain the the- 
oretical large-polaron optical conductivity spectra of 
SrTii-^NbzOa shown in Fig. 2 and Fig. 3 at 7 K and 
300 K, respectively. In each graph also the experimen- 
tal optical conductivity spectra of Ref. [17] are shown. 
It should be emphasized that in the present calculation, 
there is no fitting of material constants for the polaron 
contribution to Re<r(f2). Even the magnitude of the op- 
tical conductivity, which is often arbitrarily scaled in the 
literature, follows from first principles. 

At 7 K, the calculated optical conductivity based on 
the Frohlich model and extended for a gas of large po- 
larons as described in the present paper, shows convinc- 
ing agreement with the behavior of the experimental op- 
tical conductivity for the high energy part of the spectra, 
i.e., Ml ^ 300 meV. The experimental polaron optical 
conductivity of SrTii^Nb^Os falls down at high fre- 
quencies following the power law (derived in the present 
work and typical for large polarons) rather than as a 
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Photon energy (meV) 




Photon energy (meV) 



FIG. 2: The many-large-polaron optical conductivity compared 
x = 0.1% (a), 0.2% (6), 0.9% (c) and 2%(d). 



30 




Photon energy (meV) 



with the experiment [T7] at T — 7 K. The doping level is 



Gaussian exponent that would follow from the small- 
polaron theory. At lower photon energies Ml 200 
meV, the experiment shows distinct peaks that are not 
explained within the polaron theory. They can be due to 
other scattering mechanisms as discussed below. 

The minor deviations between theoretical and experi- 
mental Re a (Q) in the frequency range Ml ^ 300 meV 
may be attributed to the difference between the actual 
electron densities and the densities calculated on the ba- 
sis of the unit cell volume and the chemical composition. 
However, we prefer not to fit of the density. 

The optical conductivity calculated for a single large- 
polaron absorption [33] predicts an intensity 3-4 times 
larger than the experimental data for the lowest doping 
level x = 0.1%, and therefore cannot explain those data. 
For higher dopings, the overestimation of the magnitude 
of the optical conductivity within the single-polaron the- 
ory is even larger than for x = 0.1%. Therefore the 
many-polaron approach, used in the present work, is es- 
sential. 

At 300 K, in Fig. 3 (a, b ,d), the agreement between 
theory and experiment is qualitative. Both experimental 
and theoretical spectra show a maximum at the room- 
temperature optical conductivity spectra in the range 
Ml ~ 250 meV. For the doping level x = 0.9% the calcu- 
lated optical conductivity spectrum underestimates the 



experimental data, as also observed at 7 K. 



Many-body effects considerably influence the optical 
conductivity spectra of a polaron gas. First, features 
related to the emission of a plasmon together with a 
LO phonon [20] are manifested in the optical conduc- 
tivity spectra of the many-polaron gas at T = 7 K 
as separate peaks whose positions shift to higher ener- 
gies with increasing doping level. At room temperature, 
those peaks are strongly broadened and smoothened, 
and only a broad plasmon feature is apparent. Second, 
the mid-infrared optical conductivity (per particle) in 
SrTii^NbzOa is decreasing at higher doping levels due 
to the screening of the polar interactions, which is ac- 
counted for in the present approach in which S (q, t) is 
based on RPA. The effect of screening can be illustrated 
by the fact that for no ~ 10 20 cm -3 , the many-polaron 
optical conductivity per particle is reduced by about an 
order of magnitude compared to the single-polaron op- 
tical conductivity. The reduction in intensity of the po- 
laron optical conductivity band can be interpreted as a 
decrease of the overall electron-phonon coupling strength 
due to many-body effects. Correspondingly, at high dop- 
ing levels, the polaron mass m*, determined by the sum 
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rule introduced in Ref. 



7re 2 rio 



2m* 



Re (ft) dfl 



ire 2 no 
2rhi, 



(23) 



is reduced, compared to the single-polaron effective mass. 
As shown in Refs. [2011511, the sum rule [S3] remains valid 
for an interacting polaron gas. 

The large-polaron theory of the optical absorption 
based on Ref. |20j explains without any fitting parame- 
ters the main characteristics and trends of the observed 
spectra of Ref. [T7] in SrTii-^Nb^Oa, including doping- 
and temperature dependence. Nevertheless, some fea- 
tures of the experimental spectra remain to be explained. 
In particular, at T — 7 K, the pronounced peak at 
hfl ~ 130 meV in the experimental optical conductivity is 
not accounted for by the present theoretical analysis. In 
the theoretical spectra, peaks of much smaller intensity 
appear at about the same frequency. In the large-polaron 
theory, those peaks are provided by the interaction be- 
tween electrons and the LO-phonon branch with energy 
hiOL 2 ~ 58.4 meV, accompanied by the emission of a 
plasmon as described in Ref. [2D] . 

The intensity of the experimentally observed absorp- 
tion peak at Ml ~ 130 meV is considerably higher than 
described by the large-polaron theory. In the low den- 
sity limit, the experimental optical data more rapidly 



approach the single polaron limit |33| than the theoreti- 
cal predictions based on Eq. (10 1. This absorption peak 
at Ml ~ 130 meV may be provided by other mecha- 
nisms, not controlled in the present study. E. g., electron- 
phonon interaction with low- frequency non-polar (e. g., 
acoustic) phonons may contribute to the optical con- 
ductivity. The squared modulus \V q \ , which character- 
izes the coupling strength, for the deformation electron- 
phonon interaction is | V^ 2 cx q |55j . while for the Frohlich 
interaction, |V^| 2 oc q~ 2 . Consequently, for the deforma- 
tion electron-phonon interaction, the short-wavelength 
phonons may provide non- negligible contributions to the 
optical conductivity. Also, at sufficiently large q, Umk- 
lapp scattering processes with acoustic phonons can play 
a role. The treatment of contributions due to acoustic 
phonons (and other mechanisms) is the subject of the fu- 
ture work. Another possible explanation of the absorp- 
tion peak at Ml ~ 130 meV is weakened screening in the 
corresponding energy range due to dynamical-exchange 
[56]. 
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IV. CONCLUSIONS 

Many-polaron optical conductivity spectra, calculated 
(based on Ref. [20]) within the large-polaron picture 
without adjustment of material constants, explain essen- 
tial characteristics of the experimental optical conductiv- 
ity [17]. The intensities of the calculated many-polaron 
optical conductivity spectra and the intensities of the ex- 
perimental mid-infrared bands of the optical conductivity 
spectra of SrTii-^Nb^Oa (from Ref . [T7]) are comparable 
for all considered values of the doping parameter. The 
doping dependence of the intensity of the mid-infrared 
band in the theoretical large-polaron spectra is similar to 
that of the experimental data of Ref. [T7] . In the high- 
frequency range, the theoretical absorption curves de- 
scribe well the experimental data (especially at low tem- 
perature). A remarkable difference between the present 
theoretical approach and experiment is manifested on the 
low frequency side of the mid-infrared range, where the 
experimental optical conductivity shows a sharp and pro- 
nounced peak for Ml ^130 meV at 7 K. Although the 
theoretical curve also shows a feature around the same 
frequency, its intensity is clearly underestimated. This 
peak in the absorption spectrum at Ml ~ 130 meV re- 
mains to be explained. The value of the effective electron- 
phonon coupling constant obtained in the present work 
(ct e ff ~ 2) corresponds to the intermediate coupling 
strength of the large-polaron theory. 

The alternative small-polaron and mixed-polaron mod- 
els for the optical conductivity require several fitting pa- 
rameters. Furthermore, we find that the mixed-polaron 
model would need a major adjustment of the overall in- 



tensity in order to fit experimental spectra. 

Contrary to the case of the large polaron, the small- 
polaron parameters cannot be extracted from experimen- 
tal data. Moreover, the small-polaron model, for any 
realistic choice of parameters, shows a frequency depen- 
dence in the high-frequency range which is different from 
that of the experimental optical conductivity. Both the 
experimental and the theoretical large-polaron optical 
conductivity decrease as a power function at high fre- 
quencies, while the small-polaron optical conductivity 
falls down exponentially for sufficiently high Q. 

In summary, the many-body large-polaron model 
based on the Frohlich interaction accounts for the es- 
sential characteristics (except - interestingly - for the 
intensity of a prominent peak at Ml ~ 130 meV, that 
constitutes an interesting challenge for theory) of the 
experimental mid-infrared optical conductivity band in 
SrTi 1 _ a; Nb 2 ,03 without any adjustment of material pa- 
rameters. The large-polaron model gives then a convinc- 
ing interpretation of the experimentally observed mid- 
infrared band of SrTii-^Nb^Oa. 
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